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ABSTRACT 

I present the data of the shell of classical Nova Persei (1901) obtained by the 
Advanced CCD Imaging Spectrometer S3 detector on-board Chandra Observa¬ 
tory. The X-ray nebula is affected mostly by the complex interstellar medium 
around the nova and has not developed a regular shell. The X-ray nebula is lumpy 
and asymmetric with bulk of emission coming from the southwestern quadrant. 

The brightest X-ray emission is detected as an arc that covers from the west to 
the south of the central source. Part of this feature, which is co-spatial with 
the brightest non-thermal radio emission region, is found to be a source of non- 
thermal (synchrotron) X-ray emission with a power law photon index of 2.31 q g 
and a=0.68lo;)’5 at about a flux of 1.7x10 erg cm ^ s This conhrms that 
the shell is a cite of particle acceleration, mainly in the reverse shock zone. There 
are strong indications for nonlinear diffusive shock acceleration occurring in the 
forward shock/transition zone with an upper limit on the non-thermal X-ray flux 
of 1.0 xlO erg cm ^ s The total X-ray spectrum of the nebula consists 
of two prominent components of emission (other than the resolved synchrotron 
X-ray emission). The component dominant below 2 keV is most likely a non¬ 
equilibrium ionization thermal plasma of kTs=0.1-0.3 keV with an X-ray flux of 
1.6x10“^^ erg cm“^ s“^. There is also a higher temperature, kTs=0.5-2.6 keV, 
embedded, Nj^=(4.0-22.0)xl0^^ cm“^, emission component prominent above 2 
keV. The unabsorbed X-ray flux from this component is 1.5x10“^° erg cm“^ 
s“^. The X-ray emitting plasma is of solar composition except for enhancement 
in the elemental abundances (mean abundances over the remnant) of Ne/NeQ 
and N/Nq in a range 13-21 and 1-5, respectively. A distinct emission line of 
neon, He-like Ne IX, is detected which reveals a distribution of several emission 
knots/blobs and shows a cone-like structure with wings extending toward NW 
and SE at expansion velocities about 2600 km s“^ in the X-ray wavelengths. 
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The emission measures yield an average electron density in a range 0.6-11.2 
cm“^ for both of the components (hlling factor=l). The electron density increases 
to higher values ~ 300 cm“^ if the hlling factor is decreased substantially. The 
mass in the X-ray emitting nebula is (2.1-38.5) x 10“^ Mq. The X-ray luminosity 
of the forward shock ~ 4.3 x 10^^ erg s“^ indicates that it is adiabatic. The shocked 
mass, the X-ray luminosity and comparisons with other wavelengths suggest that 
the remnant has started cooling and most likely is in a Sedov phase. 

Subject headings: X-rays: stars — radiation mechanisms: thermal,non-thermal 
— supernova remnants — shock waves — binaries: close — novae, cataclysmic 
variables — stars: Individual (GK Persei) 


1. INTRODUCTION 

Nova Persei 1901 (GK Per) is one of the most extensively observed and studied classical 
nova shells over the entire electromagnetic spectrum. It is known to be a fast ONeMg nova 
with Veject — 1200 km s“^, Meject — 7 xl0“® Mq and distance ~ 470 pc (Payne-Gaposchkin 
1957; Me Laughlin 1960, Pottasch 1959). Nova Persei (1901) is the hrst recorded nova to 
show a light echo due to the reflection of light off the nearby interstellar material (Kapteyn 
1901). Goudrec (1939) shows that the light echo region and the reflection nebula detected 
later in 1917 (Barnard 1917; Oort 1951) around GK Per are results of large-scale circumstellar 
sheets of dust clouds. 

The optical remnant is 103x90 arcsec^ (diameter~0.23 pc). The images show high 
asymmetry and the remnant has evolved into series of knots and hlaments. The bulk of 
emission arises in the southwestern quadrant indicating interaction between the nova ejecta 
and the ambient gas (Slavin, O’Brien & Dunlop 1995, Seaquist et ah 1989). Furthermore, 
the remnant is detected with the Very Large Array (VLA) at 1.49 & 4.86 GHz as a non- 
thermal, polarized radio source with a spectral index of -0.67 and a flux of 24 mJy at 1 GHz 
(Reynolds & Ghevalier 1984). In addition, the flux densities are 8.7, 20.6, 29, 33 and 38 mJy 
at frequencies 4.86, 1.49, 0.608, 0.408 and 0.327 GHz, respectively (Biermann, Strom, & 
Falcke 1995). This shows the existence of shocked circumstellar or interstellar material. The 
Infrared Astronomical Satellite (IRAS) observations at 60/im and lOO/rm reveal a symmetric 
far-IR emission region extending around the nova out to 6 pc (17' either side) suggested 
to be an ancient Planetary Nebulae associated with the binary (Bode et ah 1987; Bode, 
O’Brien, & Simpson 2004). Recent optical wavelength observations of the vicinity of GK 
Per show that some of this symmetric nebulosity is produced during the quiescent mass-loss 
phase of the central binary because of the evolved nature of the secondary (Tweedy 1995). 
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More recent IR observations also indicate that the IR emission within 17' of the source is 
of material originating from the secondary (Dougherty et ah 1996). Scott, Rawlings, & 
Evans (1994) have discovered symmetric blobs of CO emission bracketing the central object 
at about ± 200 arc seconds. Overall, it is widely believed that the GK Per nebula behaves 
like a young Supernova Remnant (SNR) in the pre-Sedov phase interacting with its dense 
circumstellar medium. 

Classical novae are a subset of cataclysmic variables which are interacting binary sys¬ 
tems hosting a main-sequence secondary (sometimes a slightly evolved star) and a collapsed 
primary component, a white dwarf (Warner 1995). An outburst on the surface of the white 
dwarf as a result of a thermonuclear runaway in the accreted material causes the ejection 
of 10“^ to 10 “^ Mq of material at velocities up to several thousand kilometers per second 
(Shara 1989; Livio 1994; Starrheld 2001). Though there has been no previous detection of 
old classical nova remnants in the X-ray wavelengths (evolving like Supernova remnants), 
classical nova remnants have been detected in the hard X-rays (above 1 keV) as a result of 
wind-wind interactions in the outburst stage (Balman et. al. 1998; Mukai & Ishida 2001; 
Orio et al. 2001; Ness et al. 2003). In general, the hard X-ray radiation above 1 keV emitted 
from the hot shocked gas in nova remnants should be a powerful diagnostic for the nature 
of mass-loss mechanisms in classical nova outbursts together with the morphology, the el¬ 
emental abundances and the evolution of nova shells. A long-sought insight into this issue 
has been gained with the discovery of the remnant old shell of classical nova Persei (1901) 
in the X-ray wavelengths. 

The X-ray shell around GK Per was first discovered in a 50 ksec exposure with the 
ROSAT High Resolution Imager (HRI) (Balman & Ogelman 1999). The X-ray nebula ex¬ 
tended to 46" southwest, 60" northwest, 52" southeast, and 43" northeast of the point source 
with an elliptical shape and lumpy morphology. The count rate of the shell was found 
as O.OliO.OOl c s“^. The estimated spectral parameters of the nova shell were: Nh ^ 
1.3±0.3xl0^^ cm“^, and kT ~ 0.16T0.03 keV (~2xl0® K). The implied unabsorbed X-ray 
flux was Fa; ~ 3.0x10“^^ erg cm“^ s“^ with La, ~ 8.0x10^^ erg s“^ (470 pc). In general, what 
was detected by the ROSAT HRI were the condensations in the shell. Since the ROSAT 
HRI did not have any adequate spectral resolution, a Chandra Advanced CCD Spectrometer 
(ACIS) observation was proposed. 

In this paper, I present the characteristics of the Chandra spectrum, its components, 
and the detailed morphology of the shell. Finally, I discuss the results in comparison with the 
observations in the other wavelengths together with their implications on the novae theory 
and evolution of the classical nova shells. 
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2. THE OBSERVATION AND THE DATA 

GK Per and its vicinity is observed using the Chandra (Weisskopf, O’dell, & van Spey- 
broeck 1996) Advanced CCD Imaging Spectrometer (ACTS; Garmire et al. 2000) for a 95 
ksec on 2000 February 10, pointed 1 arc minute offset from the nominal point on S3 (the 
back-illuminated CCD) with no gratings in use yielding a moderate non-dispersive energy 
resolution. The data are obtained at the FAINT mode. Chandra has two focal-plane cam¬ 
eras and two sets of transmission gratings that can be inserted in the optical path (HETG 
and LETG; High and Low Energy Transmission Gratings). The ACIS is used either to take 
high resolution images with moderate spectral resolution or is used as a read out device for 
the transmission gratings. ACIS is comprised of two CCD arrays, a 4-chip array, ACIS-I 
(four front-illuminated CCDs); and a 6-chip array, ACIS-S (four front-illuminated and two 
back-illuminated CCDs). The spectral resolution of ACIS arrays varies depending on the 
type of CCD in use, front or back illuminated CCDs, and the row numbers in the given 
CCD as a function of distance from the read out node or the aim point. The ACIS-S3 has 
a moderate spectral resolution E/AE ~ 10-30 (falls to about 7 below 1 keV) with an un¬ 
precedented angular resolution of 0".49 per pixel (half-power diameter). The Chandra X-ray 
Center (CXC) has carried out the standard pipeline processing on the raw ACIS events, 
producing an aspect-corrected, bias-subtracted, graded (limited to grade 02346) and gain- 
calibrated event list (revision 2 data). In order to double check the standard processing, 
’’ACIS process events” thread is used to calibrate the level 1 data using the necessary gain 
and ”fef’ hies with the aid of CIAO 3.1 and the suitable CALDB (v.2.27-28). The data 
analysis is performed with CIAO version 2.1, 2.2 in the preliminary conference proceedings 
Balman (2002a,b) and in this paper with versions 2.3 and 3.1 along with XSPEC version 
11.2.0/11.3.0 and XIMAGE version 4.0 for data preparation, spectral and spatial analysis. 
Prior to the analysis, the data are also hltered from the effects of the random haring that 
occur during the observations with the ACIS-S by excluding events that are 3a above the 
mean count rate in the light curve which reduced the exposure time to about 81 ksec after 
cleaning. A preliminary analysis of the data could be found in Balman (2002a,b). 

Chandra Observatory provides good spatial resolution data with high sensitivity thus, 
an 81 ksec observation provides an excellent opportunity to detect faint X-ray sources in the 
given region of the sky. To utilize this, the CELLDETECT algorithm (Hamden et al. 1984) 
is used to detect sources in the whole Chandra held image. As a result, 16 new X-ray sources 
are detected over 6a conhdence level above the background (Kupcu-Yolda§ & Balman 2002). 
Most of the sources are found as possible Extragalactic sources. Some of them could also be 
a Galactic Cataclysmic Variable, a Low-mass X-ray Binary or a cooling Neutron Star. 
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3. THE CHANDRA IMAGE OF THE X-RAY NEBULA 

The bright X-ray nebula extends to 52" South, 41" North, 45" West, and 38" East of the 
point source with an irregular shape, most counts coming from a region centered at the SW. 
The count rate from the vicinity of the shell is O.lOiO.Ol c s“^. The count rate ratio of the 
hemisphere centered around NE to the one centered on the SW is 2:3. Figure 1 displays an 
exposure corrected image of the nebula between 0.3 and 10 keV obtained with the Chandra 
ACIS-S (S3). The ACTS transfer streak is cleaned and the pileup PSF of the central source 
is removed from the raw image. The central source (the binary system) is extracted by 
modeling a two-dimensional PSF normalized to 1.65 c s“^ (i.e., similar to the source count 
rate without pileup) using ChaRT and MARX (version 4.08) (see also section 4.1 and 4.2). 
Figure 1 shows emission from the shell that is 2a above the background with a resolution of 
0".5 per pixel. A variable smoothing is applied on the image using a Gaussian of <7=0".5-1". 
The shell is not well developed with an average surface brightness of (4.9±2.9) xl0“® photons 
cm“^ s“^ arsec“^. The X-ray nebula shows complex lumped morphology with structure (eg., 
clump/blob/hlament) sizes constrained by the smoothing performed on the image ~ 1"- 2". 
Large scale features could be collections of smaller clumps/blobs or arcs and hlaments which 
indicate radial elongations. It has a central circular region brighter on the western side and 
there is a rather faint conical emission region more extended compared with the central 
brighter part (see section 6). 


4. ANALYSIS ON THE CHANDRA SPECTRUM OF THE NOVA 

REMNANT 

Since the Chandra observation was the hrst detection of the spectrum of such a nebu¬ 
losity around a classical nova, several different models were tested on the data to understand 
the characteristics. The search showed that no single model could be fitted to the spectrum 
of the entire shell with a xl better than 8. The spectrum of the entire nova shell showed a 
thermal plasma emission below 2 keV. The harder X-ray tail of the spectrum above 2 keV 
was also best htted with a thermal plasma emission model. Equilibrium or non-equilibrium 
ionization models yielded similar hts. A power law model was ruled out for the spectrum 
of the entire remnant because of the non-physical spectral indices derived from the results 
in contrast with the non-thermal (synchrotron emission) nature of the radio data (Balman 
2002a). The spectrum was remarkably flat above 2 keV. Figure 2 shows the best htted two- 
component emission model of VMEKAL (Mewe, Gronenschild & van den Oord 1985; Liedahl 
et ah 1995) -|- PSHOGK (Borkowski et ah 1996). The shell photons were extracted using 
an annulus with an inner radius of 12".5 and an outer radius of 67".5. The background was 
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also derived from an annular region farther out around the nova normalized to the source 
extraction region. The spectrum was calculated using a minimum of 50 counts per bin and 
the channels below 0.3 keV and above 9 keV were excluded due to low statistics. 

The spectral parameters that were derived from the hts to the second component re¬ 
sulted in a very high shock temperature in excess of 50 keV and shock speed larger than 
5000 km s“^ which were inconsistent with all the measurements in the other wavelengths. 
This problem was largely due to the excess photon contribution from the wings of the PSF 
of the piled-up central source. Such a contaminating component resulted in signihcantly 
higher temperature and flux for the second component above 2 keV, and shifts in the central 
line energies (toward lower energies) of the lower temperature component below 2 keV. This 
excess emission compared with the observed PSF of the central source also needed to be 
properly cleaned in order to study the image of the shell. Sections (4.2) and (4.3) describe 
the problem and the procedure/method used to clean and remove this effect of the PSF 
wings from the shell data (i.e, images and spectra). 


4.1. The Contamination of the Shell Image by the Central Source 

Pileup results in the CCDs when two or more photons are detected as a single event. 
The consequences are spectral hardening since apparent energy is approximately the sum of 
two (or more) energies and underestimation of the true counting rate of the point sources. 
It is expected that pileup problem also distorts the PSF of the central source by creating 
holes at the centers (pulse saturation) and the PSF wings become more prominent as the 
count rates increase. As a result, appreciable central source counts can be detected in a 
region within 2' of the pointed source. This contaminates any extended emission in the 
close vicinity of the central source. The source photons of GK Per has a pileup of about 
89% where one expects the emission from the shell to be affected by the pileup PSF and 
the prominent PSF wings. In order to exclude the spectrum of the central source, GK Per, 
the nova shell photons are extracted from an annular region of 12".5 to 67".5 in radius (as 
presented in sec 4.0). This should largely avoid the contamination from the piled-up central 
source, but it will be inefficient to exclude the effects of the prominent PSF wings beyond 
12".5 out to 67".5 in radius. To segregate the counts of the nova shell from the central source 
counts in the wings of the PSF (at the vicinity of the nova shell), an archival search within 
the AGIS-S pointings (no grating or HETG observations) has been performed. The prime 
motivation of the search has been to measure the count rate in the PSF wings in an annular 
region of 12".5 out to 67".5 radius in comparison with the count rate in a circular region of 
12".5 . The particular extraction radii is assumed only to compare the GK Per observation 
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with other archival observations. The two different regions assumed are designated as: (1) 
the source dehned as the circular photon extraction region of 12".5 radius (25 pixels), and 
(2) the shell dehned as the annular photon extraction region of 12".5 to 67".5 in radius. The 
count rates within the source and shell regions are compared for several observations so that 
an acceptable estimate of the count rate in the wings of the PSF can be made for the central 
source of GK Per in the vicinity of the nova shell given the central source count rate of GK 
Per (pileup rate). Over 300 archival data of pointed sources on the AGIS-S array have been 
analyzed imposing several criteria :(1) energy spectrum of the source peaks about 1-2 keV 
and have a hard X-ray tail, (2) pointings are within a region of 2' of the nominal point on 
AGIS S-3, (3) all observations have the same frame time of 3.2 sec and obtained in the Faint 
mode, (4) no special spatial hltering has been performed during the observation for pileup 
mitigation, (5) HETG observations with bright first order data have also been excluded, (6) 
observations with pileup count rates brighter than a few c/s (i.e. 3 c/s) are also excluded 
since the scaling relation is altered, (7) observations that show strong background flares and 
anomalous high background rates have been disregarded. As a result, about 15 suitable 
pointings have been recovered that fits all the criteria. Figure 3a shows the count rate in 
a region of 12".5 radius (25 pixels), the source, versus the count rate in an annular region 
of 12".5 to 67".5 in radius, the shell, both centered at the point source. Finally, a scaling 
relation between the count rates in two different regions described as the source and shell is 
derived. The curve in Figure 3a is best htted by the quadratic function, ax^-|-bx4-c, where 
a=-0.12±0.09, b=l.ll±0.14, and c=-0.18±0.03. This has been used to approximate the 
count rate in the PSF wings of GK Per within the region covered by the classical nova shell. 
For a source count rate of 0.2 c s“^ (pileup rate), the estimated count rate in the PSF wings 
{shell region) within the classical nova shell is 0.037 c s“^, which yields an actual nova shell 
count rate around 0.063 c s“^ (instead of 0.1 c s“^). This is also checked using an archival 
AGIS-S HETG data of GK Per obtained during a dwarf nova outburst. The zeroth order 
pileup source count rate is used to estimate the count rate in the described shell with the aid 
of the derived scaling relation (same extraction radii is assumed). Gonsidering that a lower 
frame time, and spatial clipping are used to mitigate the pileup to some extend and that the 
source is in outburst, the scaling function yields similar results. I would like to note that the 
contamination problem described here could be important for Chandra AGIS observations 
of other astrophysical objects like Galaxies/Glusters and centrally hlled composite SNRs. 

In order to remove the central source from the image of the nova shell of GK Per together 
with the prominent PSF wings, a PSF is created using GhaRT. Several PSFs are merged 
to achieve the correct count rate within the nova shell region (i.e., 0.037 c s“^ in the PSF 
wings). The final central source count rate of the combined PSF in the source region is close 
to the count rate (within 8%) of the central source, GK Per, when there is no pileup (1.7-1.8 



c s“^). The un-piled central source rate is calculated using an archival ASCA data of GK 
Per (i.e., spectrum) and PIMMS. Next, the combined PSF is extracted from the total image 
to achieve Figure 1. 


4.2. The Cleaned Shell Spectra 

The central source, GK Per, has 89% pileup and has about 0.037 c s“^ in the region of 
the X-ray shell due to its prominent PSF wings. Thus, this will also affect the ’’entire nova 
shell spectrum” where the spectrum of the central source will be appearing as a component of 
the nova shell spectrum of GK Per as a result of the counts in the PSF wings. I would like to 
note that though GK Per has a dusty environment, previous X-ray observations with several 
different satellites have not recovered any dust scattering halo around the source which could 
have affected the emission below 1 keV. The result of the analysis presented in section (4.1), 
for the sources which suffer from the pileup effect on AGIS-S3, indicates that the spectra 
derived from the shell regions (as denoted in section [4.1]), basically from the PSF wings, 
may be better approximated with the spectra derived from the source regions, but scaled 
down. This is only because the spectral shapes resemble each other, and not that the spectra 
from the PSF wings have pileup photons (as in the source regions). Theoretically, pileup in 
the PSF wings should either be non-existent or very small, however the spectra from the PSF 
wings may show spectral hardening. Figure 3b shows a collection of spectra of the sources 
used for calibration where the above assumptions can be clearly investigated. The upper 
curves represent the source-region spectra and the lower curves represent the shell-region 
spectra (spectra obtained from the PSF wings). The sources have very similar pileup source 
count rates to GK Per in the stated source extraction region 0.19-0.21 c s“^ (more or less 
similar central source spectra with a hard X-ray tail - three objects are other GVs). In order 
to remove the spectrum of the PSF wings of the central source from the nova shell spectrum, 
hrst a spectrum (of GK Per) is extracted from a circular region of radius 12".5 {source-region). 
This spectrum is htted with a composite model of photoelectric absorption, bremsstrahlung 
and power law models. The best £t parameters are used to recalculate/rescale this spectrum 
using both MARX and WEBSPEG yielding a count rate of 0.037 c s“^. The parameters used 
to model this source-region spectrum are an of 1.16x10^^, a covering fraction of 0.95, a 
kT of 13.3 keV, a bremsstrahlung normalization of 5.54x10“"^, a photon index of 2.5 and a 
power law normalization of 7.93x10“^. Figure 4 shows the central source spectrum of GK 
Per derived from the source-region (upper curve), the modeled spectrum to be subtracted 
(middle curve), and the un-piled central source spectrum of GK Per (lower curve). The un¬ 
piled source spectrum is extracted from the out-of-time events of the AGIS transfer streak. 
I want to point out that all three data sets in Figure 4 are consistent with one another 
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and the un-piled source spectrum can be approximated with the modeled spectrum to an 
acceptable extend. Finally, the modeled spectrum (middle curve) is subtracted PI channel 
by PI channel from the nova shell spectrum of GK Per. Figure 5a shows the resulting total 
nova shell spectrum after removing the spectrum of the PSF wings of the central source, 
GK Per. The total number of nova shell photons between 0.3-1.6 keV is 4675 and between 
1.6-8.0 keV is 410 after removal of the background and the photons from the wings of the 
central source PSF. The shell photons are extracted using an annulus with an inner radius of 
12".5 and outer radius of 67".5 and regrouped using a signal-to-noise ratio of 13 per energy 
bin. The cleaned nova shell spectrum can not be htted with a single emission component 
with reduced of 4.0 (see Figure 5a). There is excess emission above 1.6 keV that causes 
up to J^a deviations in the residuals when the data are fitted with a single spectral model. In 
the rest of the paper, only analyses and results derived from the cleaned spectra/image will 
be presented, since there are signihcant differences between the cleaned and uncleaned data. 


4.3. The First Component (below 2 keV) 

The X-ray emission below 2 keV is found to be consistent with only the thermal plasma 
emission models. Figure 5b shows the data and the best htted two-component emission 
model of VPSHOGK-I-NEI (Borkowski et ah 1996; Hamilton et al. 1983). The shell photons 
are extracted using an annulus with an inner radius of 12".5 and outer radius of 67".5. The 
background is also derived from an annular region farther out around the nova normalized 
to the source extraction region. The spectrum is calculated using a minimum of 50 counts 
per energy bin (or a Signal-to-noise ratio of 13 per energy bin, which made no difference 
in derived parameters). The PI channels below 0.3 keV and above 8 keV are excluded due 
to low statistical quality. The spectral parameters of the nova shell are an of (0.1- 
3.7)xl0^^ cm“^, kT of 0.1-0.3 keV (~(3-4)xl0® K) and an emission measure (EM) in a 
range (0.2-54.5) x 10^^ cm“^. The noted error ranges correspond to 2a conhdence level. The 
spectral parameters derived for this component using three different composite spectra are, 
also, displayed on Table 1. The unabsorbed soft X-ray hux is Fj, (0.03-11.8) X 10-^1 erg 
cm“^ s“^ yielding an X-ray luminosity of ~ 4.3^4®2^xl0^^ erg s“^. Any two-temperature 
(TjT^Te) non-equilibrium ionization plasma model yields similar Tj and Te within 2a error 
ranges when htted to the data below 2 keV. 

The detected emission measure EM=< '^eff (after cleaning the shell spectrum) 

yields an average electron density Ue in a range 0.6-11.2 cm“^ using a volume of 4.3x10^^ 
cm^ (consistent with the X-ray photon extraction region) and a hlling factor of 1 (i.e., Ve// 
can be theoretically expressed as Ve//=1.337rv^t^f; v is the expansion velocity [see sec (1)], t 
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is the elapsed time since the eruption, and f is the volume hlling factor). If the hlling factor 
is as low as 1x10“^, then the electron density can be as high as 84-712 cm' -3 300 cm ^). 

The time scale for equipartition between electrons and ions is 


~ 2.5 X 10 


ef T, Y ( n, 

\10^K J Vl08cm-3 


-1 


( 1 ) 


(Spitzer 1978; Fransson, Lundqvist & Chevalier 1996, hereafter FLC96). The derived range 
of the X-ray temperature {Te^T^shock) and the age of the remnant implies that the emitting 
plasma is in (for ne>60 cm“^) or close to equilibrium. The ionization timescale (r=not) 
of r<1.7xl0^^ detected using the PSHOCK (non-equilibrium ionization plasma) model in¬ 
dicates it is close to equilibrium with ambient density no< 60 cm“^ (t=99 years). This is 
supported by the HI measurements in the vicinity of GK Per where the column density 
infers an ambient density of 20-30 cm“^ (Seaquist et. al. 1989). The hlling factor in the 
region where the hrst component arises is f>0.0022 calculated using uq (thus, ne=4no) and 
the range of EM derived from the spectral hts. 

The third composite model on Table 1 (i.e., §3) is constructed to derive the detected 
line energies, line sigmas and huxes together with the continuum temperature. The X-ray 
emitting plasma has mostly solar composition except for neon and nitrogen. The Ne/Ne© 
is in a range 13-21 and N/N© is in a range 1-5 (ratio of relative number fraction to H). The 
neon over-abundance is detected in a prominent emission line at around Eo=0.900-0.916 keV 
(2(7 range) corresponding to the position of the He-like Ne IX emission line triplet with a 
hux F=(1.7-11.0)xl0“^"^ erg cm“^ s“^ (three separate lines are not resolved). The range 
of the central line energy support that the (He-like) Ne IX forbidden line emission may be 
more pronounced and that the first X-ray component is an under-ionized plasma, a non¬ 
equilibrium ionization plasma, cooling toward collisional equilibrium (the nonexistence of 
the H-like neon line will strengthen this). At around 3-4x10® K, the neon is expected to 
start recombining. Nitrogen over-abundance is highly likely resulting from the He-like N VI 
emission line at Eo=0.412-0.433 keV {2a range) with a flux F=(0.35-3.16)xl0“^^ erg cm“^ 
s“^. The spectral resolution of ACIS degrades at such energies and the line prohle is smeared. 
The exclusion of the line causes deviation above three sigma level at around the line energy 
(Balman 2001). The shift in the central line energy toward lower energies is also observed 
for the N VI indicating a more prominent forbidden line emission. 

The continuum temperature derived from the thermal bremsstrahlung fit to the data 
below 2 keV is consistent with the hts performed using the other models within errors. The 
two different Gaussians used to ht the data, yield the Ne and N lines as mentioned above. 
The third Gaussian has a central energy Eo=0.549-0.562 keV {2a range) which corresponds 
to the O VH emission line. This indicates that the low temperature plasma is also rich 
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in oxygen. A fit to the data with increased oxygen abundance yields an abundance 1-3 
times the solar abundance of oxygen (ratio of relative number fraction to H). However, I 
need to caution that the O VII line overlaps with the maximum observed energy of emission 
of the low temperature component. An observation with the LETG/HETG (i.e., grating 
observations to resolve separate lines) on board Chandra to resolve the emission lines would 
be unfeasible for the shell because it has a very low surface brightness and a more recent 
HETGS observation of the central source shows no sign of the shell. 


4.4. The Second Component 

The second component is found to be consistent either with a thermal plasma emission 
or a bremsstrahlung emission in origin. The spectral parameters obtained from the fits with 
the PSHOGK model are an of (4.0-22.0) x 10^^ cm“^, a shock temperature of kT=0.5-2.6 
keV (best £t ~ 0.9 keV) and an emission measure of (1.5-160.2) x 10®"^ cm“^. Table 2 displays 
the spectral parameters for the entire shell derived from the three different composite model 
fits to the X-ray data above 2 keV. The ranges correspond to 2a confidence level and the fits 
are carried out simultaneously with the first component. Figure 5b shows a fitted composite 
spectrum. A very intriguing outcome from the spectral fits is the detection of a cold shell 
between the two components that shows high absorption from the neutral material within 
this zone. To clarify this, a fit with a single parameter is shown in Figure 5c (see also 
Table 2) where significant excess above 2 keV is observed with xl values of about 2.4 . 
Moreover, the derived plasma temperature (using thermal plasma models) from fits with a 
single absorption (Nh) parameter is very large compared with expectations. Thus, single 
absorption for both components of emission in the shell spectrum is non-physical. Fits with 
a non-thermal power law or SRGUT model also show non-physical results with positive 
spectral/photon indices (see Table 2) for a single absorption parameter for the entire shell 
spectrum. The detailed discussion on this issue can be find in section [7.1]. The ionization 
timescale (i.e., r=not) is found as r=0.9lQ^g'^xl0^'^ s cm“^. The plasma shows no resolved 
emission lines. On the other hand, the high neutral hydrogen column density absorbs the X- 
ray emission below 1.5 keV where most of the emission lines (e.g., Ne and N) would be found. 
The unabsorbed X-ray flux of the high temperature component is Fj, (0.5-25.0) x 10“^^ erg 
cm“^ s“^ consistent with an X-ray luminosity of La;=(6.6-0.2) xlO^^ erg s“^ (at 470 pc). The 
reduced values do not differentiate between equilibrium and non-equilibrium ionization 
plasma emission models. The emission measure implies a minimum average electron density 
of 2.0 cm“^ (maximum 19.3 cm“^) over the volume stated in section (4.3) with a filling 
factor f=l. If the hlling factor is decreased to about 1x10“^, the electron density increases 
to about 200 cm“^. The ionization timescale translates to ambient densities of uq < 130.6 
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cm“^ for the second component with 4.1 cm“^ calculated from the best £t results of the 
ionization timescale. Using the EM parameter values and scaling the electron density from 
the maximum limit on the ambient density (4no=ne), the hlling factor is f > 0.0002 . 

Non-thermal power law and SRCUT (Reynolds 1998; Reynolds & Keohane 1999) models 
have also been applied to the shell spectrum above 2 keV. The resulting spectral parameters 
are an of (6.0-16.0)xlO^^ cm“^ (using two different parameter), a photon index 
r=5.9, a spectral index q;= 0.32 and an unabsorbed X-ray flux F^, > 8.0x10“^^ ergs s“^ 
cm“^. Table 2 shows the detailed results of the fits using these models. The photon index 
implies steeper power law spectra not in accordance with the findings on SNRs that indicate 
existence of particle acceleration in the X-ray and radio wavelengths (see section [5.0]). Such 
values are in a range r=2-4 within errors (Rho et. ah 2002, Berezkho et ah 2003). The 
spectral index is also inconsistent with the findings on the shell of GK Per in the radio 
wavelengths (see the Introduction). 


5. THE DETECTION OF LOCALIZED NON-THERMAL X-RAY 
EMISSION FROM THE NOVA REMNANT AND SPECTRAL 

VARIATIONS 

Once the whole spectrum of the shell is determined, spectral variations on smaller scales 
are examined. The remnant is divided into four azimuthal quadrants around the central 
source from position angle PA=0° to PA=360° centered on NW, NE, SE, and SW. Figure 
6a shows a graphical representation of the extraction regions. The resulting spectra are 
simultaneously fit with a model similar to the spectrum of the entire shell. The spectral 
parameters show slight variations in different quadrants, but the variations are within 2(7 
confidence level errors in most cases. The four spectra have count rates in a range 0.019-0.006 
c s“^ (the spectrum of the PSF wings of the central source is removed with a similar method as 
described in sections 4.1 and 4.2). There are variations in each spectra at different energies 
below 1 keV (see Figure 6b). These could be related to variations in the line intensities 
(variation in elemental abundances could also be possible). NE quadrant has the least line 
intensity between 0.3-1 keV (see the 4a variation between NE and SW at around the energies 
corresponding to the O VII emission line in Figure 6b). The asymmetric morphology of the 
X-ray nebula is the result of the distribution of material in the circumbinary medium of the 
nova and geometrical effects. 

The search for spectral variations over the entire remnant shows that the spectrum of 
the nebula with two thermal components is generally consistent over the whole remnant. 
On the other hand, the radio synchrotron emission from the nova shell strongly suggests 
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that non-thermal emission associated with the accelerated electrons should be recovered. 
Such a component that is not prominent enough to be detected from the entire shell, then 
could be localized within the shell. Most of the efforts to derive spectral variation across the 
remnant suffered from the low surface brightness of the shell except for the highest intensity 
zone at the edge of the remnant which is an arc of emission covering from NW to S of the 
X-ray nebula. This is also co-spatial with the brightest radio synchrotron emission region 
at 1.4 GHz (see section [6.1]). The spectrum of the brightest X-ray emission region, Peak 
Spectrum, is extracted from the edge of the shell using a sector that subtends about 70° to 
the central source with a width of nine arc seconds (the region denoted as (2) in Figure 6a). 
Next, the Peak Spectrum is cleaned using a similar procedure described in section (4.1 & 
4.2). In addition, the rescaled source spectrum is corrected further using the total number 
of central source photons within the extraction annulus of the entire nova shell spectrum (~ 
3000 counts) and a radial distribution of r“^ (i.e., N(r)~Ar“^ see Chandra calibration pages 
for further details). The simulated central source spectrum (as described in section 4.2) is 
subtracted PI channel by PI channel from the Peak Spectrum. The cleaned Peak Spectrum 
is regrouped with a minimum of 20 counts per energy bin. The uncleaned Peak Spectrum 
has a count rate of 0.0099±0.00072 which reduces to 0.00925±0.00057 after the cleaning 
process is completed. It is apparent that the cleaning process is not crucial and necessary 
for the bright X-ray emission region since this region is far away from the central source and 
small enough so that the contribution from the PSF wings of the central source is marginal 
(below the local background). The total number of photons between 0.3-1.6 keV is 670 and 
between 1.6-8.0 keV is 105 after the removal of the background and the photons resulting 
from the wings of the central source PSF. Because of the low statistical quality PI channels 
below 0.3 keV and above 8.0 keV have been excluded. 

An investigation of the cleaned Peak Spectrum, performing hts to the data reveals that 
it has two components. Table 3 displays all the results from hts with several composite 
thermal and non-thermal emission models. Figure 7a shows the Peak Spectrum htted with 
a single model of emission where the reduced xt value is 3.4 and the excess in the band 
1.6-8.0 keV is signihcant. Thus, existence of two components of emission is justihed. The 
hrst component resembles to that of the hrst component of the entire shell with similar 
spectral parameters. The spectrum of the second component from the X-ray Peak region, 
this time, yields consistent hts with the SRCUT or power law emission models as opposed 
to a thermal plasma model which yields non-physical results. All the hts thermal plasma 
models result in temperatures in excess of 7 keV (2(T lowest limit and best ht result ~ 20 
keV) which is inconsistent with the expectations from this remnant. The hts also reveal 
absorption diherence between the two components in question and the hts involving a single 
absorption parameter for the two components yield signihcant excess in the 1.6-8.0 keV 
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range with reduced xl value of 3.0 . Figure 7b and Figure 7c display examples of the htted 
Peak Spectrum that has a two-component model with two different parameter, and a 
two-component model with a single parameter for both components, respectively. The 
spectral parameters for the second component derived using a simple power law model are 
an Nh of cm“^, a photon index of F=2.3lQg and an unabsorbed X-ray flux 

of 1 . 71 ^ 4 '^ X10“^^ 6^8 cm“^. The errors correspond to 2(j confidence level. The X-ray 

flux of this non-thermal component translates to an unabsorbed X-ray luminosity of about 
4.6x10^° erg s“^. More appropriate models for synchrotron emission (from SNRs) involve 
exponentially cut-off power law distribution of electrons in a magnetic held (Reynolds 1998; 
Reynolds & Keohane 1999). This is largely because the synchrotron losses or escaping 
particles are expected to decrease the hux density at the X-ray energies falling below the 
extrapolated radio spectrum. The cleaned spectrum from the X-ray Peak region (region [2] 
in Figure 6a) is htted with the SRCUT model for a more proper treatment. The resulting 
parameters are : N/f=5.0^4;7Xl0^^ cm“^; q;=0.681o; 45; and z/breaA:=l-2lo!3xlO^* Hz. I would 
like to make a note that since the photon extraction region for the Peak Spectrum is very 
small and farther away from the central source, the cleaning process have slightly changed 
(i.e., corrected) the best ht parameters leaving the error ranges similar/unaltered (for both 
SRCUT and the power law models). The normalization parameter, which is the radio energy 
hux at 1 GHz (in Jy), is hxed at the calculated value of 0.024 Jy (Reynolds & Chevalier 
1984; Biermann et al. 1995). The energy hux and the unabsorbed X-ray luminosity is the 
same as the values calculated for the simple power law model. The ht successfully recovers 
the correct spectral index ,«, for the radio spectrum (q;= 0.67; see Introduction). This is also 
consistent with the range of X-ray spectral photon index. Thus, emission from the same 
electron population has been detected as the extension from the radio wavelengths. 

It is expected to detect gamma-ray emission as well as X-rays from the shell-like SNRs 
primarily due to bremsstrahlung of accelerated relativistic electrons oh pre-shock material 
together with contribution from an Inverse Compton emission (oh of Cosmic Microwave 
Background), the pion decay interactions, and the secondary p-p interactions (Berezhko & 
Volk 1997; Grassier et al. 1998; Stunner et al. 1997; Baring et al. 1999, Berezkho et al. 
2003 and references therein). The critical frequency (also the break frequency in the above 
paragraph) for synchrotron emission from an electron of energy E is 


^break 3.2 


Bsind 
2 X lO-^G 



2 

GHz 


( 2 ) 


(Chevalier 1999). Using the cut-oh (break) frequency derived from the SRCUT model hts, 
I estimate E^aa; — 15-30 TeV for the maximum energy of accelerated electrons radiating at 
(0.4-1)xlO^® Hz for the remnant of GK per (B=75-47/iG; Seaquist et al. 1989). 
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A spectrum is extracted from the outer ridge (outer ridge: region denoted as (3) in Figure 
6 a) using a PIE extraction region with a width of nine arc seconds (the same size annulus- 
sector as the photon extraction region of the Peak Spectrum). A £t with the SRCUT model 
of emission yields an Nj:/=2.2]'^q;®x 10^^ cm“^, Q;=0.76]'^Q;g2 and z/fereaA:> 1-8x10^® Hz (x^=1.17 
(37)). The unabsorbed X-ray flux from the outer ridge is erg s“^ cm“^. The 

spectrum of the inner ridge (inner ridge: region denoted as (1) in Figure 6 a) is extracted, 
using a similar photon extraction region as in the derivation of the Peak Spectrum and the 
outer ridge (with a width of nine arc seconds). A fit with the SRCUT model of emission 
yields an Ni^= 2 . 1 ]'^^; 3 XlO^^ cm“^, q;= 0 . 581 o;o 7 and Hz {xl=0.7 (23)). 

The nnabsorbed X-ray flux from the inner ridge is l.llg 07 X 16 “^^ erg s“^ cm“^. All the 
spectra are cleaned (for the wings of the PSF of the central source) using a similar procedure 
to the Peak Spectrum where the amount of cleaned PSF photons is largest for region (1) and 
lowest for region (3). The fits with the SRCUT model for the regions 1, 2 and 3 in Figure 6 a 
yield spectral indices consistent with the radio observations of the remnant which are fonnd 
to be in a range 0.55-0.95 (between 1.4-4.9 GHz; see Seaqnist et al. 1989). The radio spectral 
index is variable over the remnant where it is steepest at the outer ridge and more flat in the 
interior regions. The normalization of the SRCUT model fits is fixed at the calculated value 
of 0.024 Jy at 1 GHz as mentioned earlier. Thus, it is clear that there are three components 
in the spectrum of the X-ray nebula of GK Per: A cool thermal component associated with 
the forward shock/transition zone that shows line emission; a deeply embedded hot thermal 
component without emission lines and a ’’localized” embedded non-thermal power law (i.e.. 
Synchrotron) emission component coming from the region where the shell interacts with the 
circnmbinary medinm aronnd the nova. An elaboration of this scenario is in section [7.1]. 


6. THE Ne IX EMISSION LINE IMAGE AND DETAILED IMAGING OF 

THE REMNANT 

The nebula displays a different structnre compared with the broad band Chandra image 
at the emission line energy detected aronnd 0.907 keV. The Ne IX emission line image shows 
a symmetric cone-like shape centered at the SW direction (see Figure 8 ) and the flattening on 
the SW is more prominent. The image is prepared by extracting the energy channels (0.8-1.0 
keV) forming the full width of the line at the continnum level and snbtracting on pixel-to- 
pixel basis a suitable continuum image from the line image using PI channels corresponding 
to 1.12-1.4 keV. This continuum image, also, subtracts out any central source photons since 
it includes both the source and the shell continnum in the channel range. I need to note 
that some of the line emission and diffuse neon emission could have been lost in this process, 
bnt Figure 8 is a cleaned image. Both of the line and continnum images are regrouped by 



two pixels and smoothed by a Gaussian of <7=2" before subtraction. The resolution of the 
hnal image (continuum subtracted) is 1" per pixel and shows emission from the shell that is 
2a above the background. The neon image of the shell clearly reveals a cone, a one sided 
shell (i.e., hemisphere), where the remnant rams into a wall like structure and is compressed 
in that direction. The image also shows a very lumpy morphology with clump sizes limited 
by the resolution of the image ~ 2". The size of the cut-cone is ~ 63" in the flattened 
zone and ~ 112" -116" in the wings toward the North and East in the image. Assuming 99 
years for the elapsed time and a linear expansion law, the velocities required to create the 
wings are ~ 2600-2800 km s“^ and the flattened zone is ~ 1100 km s“^. The high velocity 
necessary for creation of the ’’wings” detected in the Ne IX image are more than twice the 
expansion speed of the nova shell detected in the other wavelengths and about the same in 
the SW region which is the flattened zone. This shows that the ejecta material is mixed 
in the forward shock. Using the extensions above and the law of cosines, the cone angle is 
calculated to be ~ 74°. The axis of the cone is in the same direction with the axis of the 
polar cones attributed from the material lying in the light-echo region (Seaquist et al. 1989). 

The reverse shock regions are prone to formation of knots because of Rayleigh-Taylor 
instabilities. The width of the interaction region between the forward and reverse shocks 
shrinks considerably with increasing compression ratio (due to existence of nonlinear diffusive 
particle acceleration) and the convective instabilities can reach all the way into the forward 
shock. As a result, clumps and hlaments can be found in the vicinity of the forward shock 
region. For low compression ratios, any enhancement of the radial magnetic held could play 
the same role and Rayleigh-Taylor instabilities can still extend all the way into the shock 
front (see also Blondin & Ellison 2001 for a discussion). This explains the existence of knots 
of emission (neon knots) in the forward shock zone of GK Per and also mixing of the ejecta 
material into the forward shock. However, this does not necessarily account for the large 
asymmetry detected in the remnant. 

Therefore, a contributing factor to the above argument for the structure in the emission 
line would be the existence of undecelerated, cool and dense neon knots in the shell that 
have caught up with the forward shock and surpassed it in SE and NW directions where 
the interstellar density is low. Such a phenomenon was also suggested for the knots in Gas 
A, Vela, Gygnus Loop and others. The SE-NW direction is also the site for the reflection 
nebula, the symmetric far-IR emission nebulosity, and suggested to be a region of low density 
(Lawrence et al. 1995; Tweedy 1995; Bode et al. 1987; Seaquist et al. 1989). This is 
also supported by the GO map of the vicinity of the nova which shows a relatively empty 
environment toward the north and east of the nova and a prominent GO emission toward 
the south (Hessman 1989; Scott et al. 1994). Moreover, detected neon clumps could be 
formed at an early stage in the eruption of the nova and matter could be ejected as knots in 
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high velocity winds during the early outburst stage (Shore et al. 1997). This also implies a 
binary orientation origin for the asymmetry in the remnant together with the non-uniform 
density prohle of the circumstellar medium around the nova. 

In addition to the Ne IX emission line image, other possible emission line images were 
constructed at around 0.58 keV (O VII line) and 0.43 keV (N VI line). The resulting images 
indicate that the shape of the nebula at these energies do not indicate a special structure 
and the wings seen in the neon image are not detected. Further energy dependence within 
the X-ray nebula is investigated by constructing images at different energy bands. Figure 9 
displays the image of the X-ray nebula in three different energy ranges: (a) 0.25-0.50 keV, 
(b) 0.51-1.50 keV, and (c) 1.51-8.0 keV. The hrst two images show emission from the X-ray 
nebula that is 2(T above the background at 1" per pixel resolution, smoothed using a Gaussian 
with a variable width a=l"-2". They are also corrected for the contamination caused by the 
ACIS transfer streak and the PSF wings as described in sec. (4.1). All three bands indicate 
different structure. The first two bands (up to 1.5 keV; Figure 9a,b) show generally a hlled 
semi-circular structure brighter on the western side with the brightest zone at the flattened 
rim which is more evident in the second band. This resembles the structure of the total band 
image (see Figure 1) except for the ’’wing” structure detected predominantly in the neon 
emission line (seen also in the second band image). The second band image is more cone-like 
owing largely to the motion of the neon knots as described in the previous paragraphs. The 
third band (hardest X-rays) image is partially resolved (detected spectrally in sec. 4). The 
original discovery of the X-ray nebula of GK Per was also at a resolution of low signihcance 
with the ROSAT HRI between 0.2-2.4 keV (see Figure 1 in Balman & Ogelman 1999). 
The Ghandra observation has recovered most of the emission and structure missed out by 
the ROSAT HRI, however, the emission above 2 keV is still not completely resolved, but 
resembles the structure of the ROSAT HRI detection and, also. Figure 9a. Figure 9d is 
constructed to stress the resolved emission regions at and above 3a conhdence level with a 
lower resolution of 3" per pixel. It is smoothed using a Gaussian with a width of a=Q''. The 
most prominent region revealed in Figure 9d in the 1.51-8.0 keV band is ellipsoidal with a 
major axis of 55" and a minor axis of 27" coincident with the brightest region of the shell 
in 0.2-1.5 keV range as shown by the overlay of contours in Figure 9d. Part of this elliptical 
emission zone is also the region where the Peak Spectrum is derived and the non-thermal X- 
ray emission is detected. The different structures in the three different energy bands support 
that the X-ray emission from the shell of GK Per has different components. 



6.1. Comparisons with other Wavelengths 


In order to investigate the conditions within the shell and determine the evolutionary 
stage of the remnant, spectral and spatial comparisons in different wavelengths are necessary. 
Detailed results from the recent data on GK Per are not available. However, a preliminary 
comparison of the X-ray image has been made by aligning the image with the recent images 
obtained in the optical and radio wavelengths. Figure 10a shows the X-ray image and an 
overlay of contours derived from the HST image obtained with the [Nil] hlter (brightest 
emission line in the optical wavelengths) (HST image: Shara 2002, private communication). 
The combined image reveals that the [Nil] knots and hlaments are in general coincident with 
the central circular region of the X-ray nebula, not with the ’’wings” detected predominantly 
in the Ne IX image (see also Figure 9). Both images do not show a well developed shell and 
are brighter on the Western hemisphere rather than the East. This is not an effect of the 
X-ray absorption since the spectral hts do not yield signihcant variation in the Nh parameter 
in different directions within the nebula. In general, the [Nil] image is larger toward the SW 
of the nova which indicates cooling in the forward shock and that the remnant has evolved 
in comparison with the aligned images in Seaquist et ah (1989). The optical shell, mostly 
composed of knots and hlaments, suggests large-scale density gradients in the vicinity of 
the classical nova. Though the contours look circular in geometry, the HST data indicates 
an elliptical velocity gradient in the shell with the fastest material lying in the NW to SE 
direction (Shara 2002, private communication), which is the direction the elongation in the 
X-ray nebula lies, as detected in the image of the Ne IX line and Figure 9b. 

Since there is X-ray emission both from the reverse and forward shocks, the line emission 
in the UV and optical wavelengths coming from within the shell and the pre-shock zone 
should be affected by the ionization due to X-ray irradiation (i.e., particularly, the luminosity 
and evolution of the Ha emission line; Chevalier & Fransson 1994; Chevalier & Oishi 2003). 
The [OHI] emission line is a characteristics of the shock zones in circumstellar interaction, 
however the [OHI] emission from the shell is weak and mostly suppressed in the SW region 
where the X-ray emission is the brightest (Seaquist et al. 1989; Slavin et ah 1995). This 
could be because the expanding shell is stopped very efficiently and the shock speeds are low 
in this region (shifting the emission to Nil and Ha lines). The densities of the optical knots 
in the SW are also detected to be high in comparison with the ones in the NE, SE and NW, 
thus the prominent emission lines are that of [Nil] and Ha (Seaquist et al. 1989). 

The image of the X-ray nebula has been aligned also with a recent VLA image obtained 
at 1.425 GHz in 1997 (VLA image: Seaquist 2002, private communication). Figure 10b 
shows the X-ray image of the shell with an overlay of radio contours. The radio image is 
mostly circular in structure coinciding with the central circular part of the X-ray nebula 



19 


and is slightly smaller than the X-ray image (as compared at a level of 2a emission, see also 
Balman 2002a). The flattening on the SW is evident in both images. An intrigning resnlt of 
this combination image is that the peak of the non-thermal radio emission falls slightly offset 
toward inside of the shell from the brightest X-ray region (X-ray peak of the nova shell). 
This indicates that the particle acceleration conld be occnrring at the reverse shock zone 
consistent with the high absorption associated with the non-thermal X-ray component (see 
section 5.0). In addition, the onter ridge of the radio shell falls in the X-ray bright zone, as 
well. An example similar to this has been detected in the Chandra observation of the yonng 
snpernova remnant IE 0102.2-7219 (Gaetz et al. 2000). Moreover, in Fignre 10b the radio 
wavelength contonrs show a trace of the ’’wings” at a 1.5a conhdence level (Fignre 8, see also 
the second band image in Figure 9). The part of the radio wavelength (1.425 GHz) image 
associated with the ’’wings” seems to be a part of a much larger and fainter emission which 
may not be associated directly with the shell, but the circumstellar/interstellar environment 
of the nova. It extends from NW to SE which is in accordance with the lower circumstellar 
density region and the location of the PN and/or the reflection nebulosity mentioned earlier. 
The content of the radio spectrum along the ’’neon wings” is at this time unknown. 


7. DISCUSSION 

7.1. On the Origin of the X-ray Nebula 

In general, the Chandra observation of the shell of nova Persei (1901) has revealed that 
the remnant evolves similar to the young supernova remnants expanding into its circumstellar 
medium, a leftover wind ejecta expelled in the course of the binary evolution, and a GO 
enhanced interstellar medium in the vicinity (similar to Type II SNRs: Ghevalier & Fransson 
1994 hereafter GF94; see also Truelove and Mckee 1999 for a general review). A conspicuous 
interpretation of the spectrum of the nova shell is such that the origin of the emission below 
2 keV is partly of the ejecta (e.g., Ne and N emission lines) and the circumstellar matter, 
namely the forward shock (also the transition zone), whereas the emission above 2 keV is of 
the reverse shock region (owing to the absorption by high equivalent neutral column density). 

There is signihcant difference (Sa) between the neutral Hydrogen column densities of the 
two components of the nova shell spectrum. This could be explained in the context of SNR 
evolution where such a cool layer may exist between the forward and the reverse shocks in 
the early phases. The size (Ar) of this cold shell is less than 5x10^^ cm (less than 0.0005 arc 
seconds at 470 pc distance) for an SNR 100 years old (GF94), with a range of density 2x10® 
to 3x10^^ cm“® (GF94). Multiplying the density range with the possible size of this region 
yields a value of Nh in the cold shell in a range Ixl0^^-2xl0^^ cm“^. The column density 
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in the cool layer decreases with time and can be expressed as Ncooz—9x10^^(n-4)(t/l 1.57 
days)“^ cm“^ (FLC96). I have assumed 1200 km s“^ for the expansion speed, a power law 
index of 2 for the circumstellar density gradient, an M=1x10“^Mq yr“^ for the circumstellar 
wind and a wind velocity of 10 km s“^. For a range of n from 10 to 200 (n=power law index 
for the density gradient in the ejecta), 99 years of age yields Niif~1.2xl0^‘^-8.7xl0^^ cm“^ 
in the cool layer. SNRs at the age of 100 yrs are about 2 pc in radius (4 pc in size) whereas 
the radius of GK per is only 0.11 pc. The size of an SNR at the same evolutionary age 
with GK Per is larger by a factor of 20 at the least. Thus, time may disperse a cool shell 
between the forward and reverse shocks in SNRs of high velocities and large sizes, but, since 
GK Per is considerably slower (compared with an SNR 100 years old) and more compact 
with respect to an SNR at the same evolutionary age, this cool absorbing layer could stay 
intact. Moreover, the size of GK Per is equivalent to an SNR a few years old. In such cases, 
a distinct cool layer is likely (if the reverse shock is radiative) and the emission from the 
reverse shock would be absorbed. 

The high X-ray absorption detected for the shell of GK Per between 0.5-2 keV is not 
largely due to the neutral hydrogen, but to the neutral metals owing to the fact that the 
absorbing layer is the cooled part of the reverse shock/ejecta itself. The composition of 
the ejecta is enhanced in metal abundances of N, Ne and possibly O. Even the circum¬ 
stellar wind material is at least of solar composition. The absorptivity per hydrogen atom 
greatly increases above 0.1 keV due to the existence of metals like N, O, Ne, and Fe of 
solar abundances (Wilms, Allen & McGray 2000 and references therein) which would yield 
large equivalent column densities of neutral hydrogen. The photoelectric absorption cross- 
section ax oc (Z=atomic number) and the optical depth increases since r oc In 

order to check this, the same composite models in the Tables 1, 2 and 3 are htted using the 
multiplicative models VARABS and VPHABS (Balucinska-Ghurch & McGammon 1992) in 
XSPEG to model the absorption in the second component. These models allow calculation 
of equivalent column densities of neutral Hydrogen together with elemental abundances of 
the absorbing medium. The results show that the X-ray data is consistent with a range of 
equivalent neutral Hydrogen column density of cm“^ attained by, in general, 

a plasma of solar abundances with N/N 0 =lO. 9 l 9 ^^'^ Ne/Ne 0 = 17.9lgQ. The elemental 

abundances determined from the cool layer are slightly higher (esp. N) than what is recov¬ 
ered from the analysis of the first component which is the forward shock mixed with the 
ejecta. Assuming that the cool layer portrays a better representation of the ejecta, then the 
elemental abundances of Ne and N could be higher than displayed on Tables 1 and 3. The 
total column density of metals causing the absorption in the cool shell is estimated as N^< 
2 x 10 ^° cm“^ which yields a total mass in the cool layer of McooZs/ieZZ<l xlO“®M 0 (less than 
1.5% of the ejecta mass). 
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The reverse shock speed (using the temperature derived for the second component) 
is found as 872 km s“^ and the temperature of the hrst component (forward shock) im¬ 
plies shock velocities of 400 km s“^ for the forward shock calculated using the general rela¬ 
tion kTg = , assuming Rankine-Hugoniot jump conditions in the absence of 

cosmic-ray acceleration. The average expansion speed detected for the nova is about 1200 
km s“^ in the optical and UV wavelengths (see also sec 1). This value is consistent with the 
values obtained by measuring the size of the X-ray nebula in the SW direction (see section 
6.0), but inconsistent with the speeds calculated for NW and SE directions (see section 6.0). 
The shock speed calculated using Rankine-Hugoniot jump conditions for the hrst component 
associated with the forward shock is also not in accordance with the expansion speed of the 
shell, at all. 

The reverse and forward (circumstellar) shock temperatures scale with Tj-ev = ((3 — 
sY/{n — 3)^) Tcs, where s is the power law index of the circumstellar density gradient and 
n is the index for the radial dependence in the ejected material as noted earlier (FLC96). 
Following this, the ratio of the forward and reverse shock temperatures for the nova remnant 
are about 4-16 and the shock velocities are about 2-4 (i.e. s=l-2 and n=7). As discussed in 
the previous paragraph, n is larger than 8 (for a radiative reverse shock) and thus the actual 
contrast between the forward and reverse shock temperatures are expected to be larger than 
the above value. This approximation yields a minimum estimation of about 1600-3200 km 
s“^ for the forward shock speed (~ 8 keV) using the computed reverse shock velocity. There 
is no component detected in the X-ray spectrum with the stated shock speed range above and 
the corresponding temperature. Therefore, the X-ray temperatures does not follow from the 
shock speeds/expansion speeds as expected from simple shock evolution. Some of the detected 
electron temperatures for the blast wave interaction region in young SNRs are found to be 
as low as ~ 10®“'^ K (e.g., IE 0102.2-7219: Hughes et ah 2000, Ellison, Slane & Gaensler 
2001; Kepler: Decourchelle et ah 2000; SN 1006: Bamba et ah 2003). This range of 
shock temperatures are at least a factor of 10 less than what would be expected from such 
young SNRs. Recent theoretical calculations show that the existence of nonlinear diffusive 
shock acceleration (particle acceleration at the shock zones) modihes shocks, changing their 
structure and evolution (Decourchelle et al. 2000; Ellison, Decourchelle & Ballet 2004 and 
references therein). The particle acceleration at the shock zones (forward or reverse shocks) 
causes the Mach number to decrease in the post-shock (i.e., sub-shock) regions making the 
shocked gas more compressible (increasing the compression ratios) resulting in a transition 
region between the reverse shock and the forward shock that is narrower and denser. As 
the efficiency of the acceleration increases, the post-shock temperature (esp. for protons) 
cools significantly (e.g., factor of 10; Ellison et al. 2004 and references therein). Therefore, 
existence of particle acceleration modihes the shocks and Rankine-Hugoniot conditions are 
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no longer valid. Consequently, the low temperature (0.1-0.3 keV) and the inconsistent shock 
speed (400 km s“^) calculated for the hrst component of the entire shell spectrum associated 
with the forward shock could be well explained by nonlinear diffusive shock acceleration. It 
is important to note that the higher densities and lower temperatures will also lead to more 
rapid Coulomb equilibration than expected in cases without acceleration. 

Considering both of the components and derived electron densities in sections (4.3) and 
(4.4), the mass in the X-ray emitting nebula is estimated as (2.1-38.5) x 10“^ Mq using the 
simple relation, M^heii — '^eff- The calculated shocked mass is consistent with the 

ejecta mass expected in classical nova explosions (see sec[l]) and the range is 10-50 times 
the ejecta mass. The expansion speed for nova Persei 1901 yields a total kinetic energy of 
(1.1-6.2) xlO^^ ergs for the ejecta (1/2 Mg^vT and for M^j and vT see the Introduction). The 
total explosion energy of a classical nova event is of the order of 10^®-10^^ ergs (i.e., binding 
energy of the envelope of the white dwarf; Eo= 2GyiwD^envl^wD)- Assuming that the shell 
will be radiative at the time maximum X-ray luminosity is reached, then Lmaa;—Eo/tmaa;- If 
tmaa; Were 99 years (time since the eruption), we can calculate that Lmax—10^^“^® ergs s“^. 
The kinetic energy dumped in the ejecta (about 1% of the total explosion energy) translates 
to about Lmaa:—10^^ ergs s“^ for a radiative shell. The reverse shock luminosity is expected to 
decrease in time by a factor of 100-1000 times after 100 yrs (CF94) and the luminosity should 
be, then, in a range ergs s“^ which is in accordance with the detected luminosity of 

the second component (the reverse shock; La;=(6.6-0.14)xl0^^“^^ ergs s“^). On the other 
hand, given the low luminosity of the shell originating in the forward shock a few x 10^^ 
ergs s“^ and the fact that the calculated maximum shocked mass is larger than the ejecta 
mass, the remnant is most likely evolving in the Sedov phase as an adiabatic remnant. 

It is of evolutionary importance to note that the environment of GK Per resembles the 
class of SNRs that interact with molecular clouds (Seaquist at al. 1989). The Figure 1 and 9 
reveals that the shell has a lumpy structure which could, also, support the existence of cloud 
crushing; shocks with the molecular clumps (White & Long 1991). A part of the emission of 
the second component could be due to this. Interactions of SNRs with dense environments 
have been modeled and old remnants like W44 and IC443 are recovered as the candidates 
for such SNR evolution with a colder radiative shell in the X-rays and a hot isothermal gas 
around 1 keV at the center (Rho & Petre 1998; Shelton et al. 1999; Chevalier 1999; Bykov 
et al. 2000). The modeling of SNR W44 (interacting with a dense inhomogeneous region 
in the vicinity of a molecular cloud) yields a one sided X-ray emitting hemisphere, a hot 
interior and a colder radiative shell as modeled by Cox et al. (1999). This has a remarkable 
resemblance in structure to the remnant of GK Per except that the outer shell is not radiative 
and cold due to age, but rather adiabatic and cold most likely as an effect of the particle 
acceleration and modihed shock evolution. Overall, the nova remnant of GK Per could be a 
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younger remnant that will eventually resemble older remnants like W44 or IC443. 

The photon emission in the X-ray wavelengths presumably arising from the acceler¬ 
ated particles could be a non-thermal bremsstrahlung process resulting from a particular 
non-Maxwellian distribution of accelerated electrons (Asvarov et ah 1990) or synchrotron 
emission from a power law distribution of such electrons. The non-thermal synchrotron 
emission from the nova remnant as reported in section (5.0) is coming from the brightest 
X-ray emission region which is most likely an X-ray filament. Lately, such non-thermal £1- 
aments have been recovered in SNRs with similar photon indexes in a range 2.1-2.5 (e.g., 
Bamba et al. 2003; Uchiyama et al. 2003, Ueno et al. 2003). The spectrum of accelerated 
particles in the X-ray band should be highly sensitive to the local neutral density and have 
a non-thermal/thermal bremsstrahlung component at high no for hard X-ray emission. The 
non-thermal bremsstrahlung intensity is 1.0xlO^®nNoE“’"/(r-l) cm“^ erg“^ s“^ (cf. Cheva¬ 
lier 1999). The emission from radio to the X-ray band is dominated by synchrotron emission 
from energetic electrons for B > 10 fiG with bremsstrahlung dominating X-rays for B < 3/iG 
(see discussion in Ellison, Berezhko, & Baring 2000 and Ellison et al. 2004 for more details). 
The estimations on the circumstellar magnetic field intensities in the non-thermal radio peak 
region and the following outer ridge of GK Per are 75 and 47 fiG, respectively (Seaquist et 
al. 1989). As a result, one expects a synchrotron X-ray emitting region co-spatial with the 
peak and outer ridge detected in the radio wavelengths of GK per as presented in this paper. 

The X-ray luminosity and the power law index from the brightest X-ray (peak) region 
can be used to determine the non-thermal electron density and the total energy in the 
accelerated electrons. The accelerated electron population will have a distribution in basic 
terms Ne(E)=K E“^ with an individual electron losing energy via synchrotron process as 



averaging over all pitch angles and taking electron velocity ~ c, where ax is the Thomson 
cross-section. The electron distribution described above will have a multiplicative exponen¬ 
tial decay term to explain the age-limited, and necessarily the loss-limited or escape-limited 
SNRs more properly (see Reynolds 1998 and references therein). Following from the pre¬ 
scriptions of Reynolds (1998), the ratio of the maximum electron energy for the remnant 
of GK Per and the theoretical energy (i.e., Ej in eq. 27 and 26 of Reynolds (1998)) of an 
initially infinitely energetic electron in the same magnetic field (as in the remnant of GK 
Per) radiating for the same duration time t predicts that the remnant is age-limited (young 
for synchrotron losses to dominate or particle escape to be prominent). 

The injection efficiency for the remnant {rf = Unon-thermai/nthermai) is proportional to the 
ratio of the electron number densities of the non-thermal and thermal electrons. I estimate 
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that 7 ] = 1.4x10“^ for the brightest X-ray (peak) region of the remnant of GK Per. This is 
about the same order of magnitude with the r/ reported for SN 1006 (Bamba et al. 2003) 
which is the most efficient SNR known to accelerate particles to very high energies close to 
the knee. The non-thermal electron density is calculated as T^non-thermai = 3.01x10“^ cm“^ 
for a hlling factor of 1 {nthermai=‘^-0 cm“^ for f=l). A basic power law electron distribution 
is assumed for the calculation, along with the integration limits : (1) a minimum energy (the 
injection energy) of 1 keV (the temperature of the thermal plasma of the second component), 
(2) a maximum energy of 15 TeV. In accordance with the photon extraction region, a volume 
of 2.8 xlO®^ cm ^ is also, assumed. 

The non-thermal energy density of the brightest X-ray emission (peak) region is es¬ 
timated as 3.0x10“^*^ erg cm“^. This is calculated using the same integration limits and 
volume described in the above paragraph together with a distribution of power emitted by 
the given power law electron distribution (i.e., [dP(E)/dE]~i?[dN(E)/dE]). The correspond¬ 
ing energy densities of the magnetic field and the thermal plasma are 2.3x10“^*^ erg cm“^ 
(uB=B^/87r ) and 5.1x10“® erg cm“^ (utftermaZ=3/2 UekT), respectively. The energy densi¬ 
ties are a factor of 10 higher compared with SN 1006 by Bamba et al. (2003). The total 
energy in the thermal plasma, the magnetic field and the non-thermal electrons are then, 
2.2x10^^ ergs, 9.5x10^® ergs (if total volume is assumed and 6.2x10^^ ergs if the photon 
extraction region is assumed), and 1.2x10^^ ergs (if total volume is assumed and 8.4 xlO*^^ 
ergs if the photon extraction region is assumed), respectively. As a result, the magnetic held 
is almost/in equipartition with the non-thermal electrons (taking that these calculations 
are estimations). The non-thermal protons should carry more energy than the non-thermal 
electrons (in a range 1-2000 times the total energy of electrons [2000~mp/me]). 

There are strong indications of possible nonlinear diffusive particle acceleration occurring 
in the forward shock zone as well as the reverse shock zone; (1) the unexpectedly cool X- 
ray temperature, (2) the strongly compressed interaction region between the forward and 
reverse shocks where the ejecta material is detected to be mixed in. However, including a 
second power law or SRCUT model of emission in the hts applied on the X-ray data below 
2 keV does not yield physical or statistically consistent results unless the hrst component 
is htted with a bremsstrahlung model and three Gaussians for emission lines along with 
the SRGUT/power law model. Inclusion or exclusion of the SRGUT/power law model does 
not affect the reduced of fh® fits. Thus, an upper limit on the non-thermal X-ray flux 
from the forward shock is Fa;<1.0xl0“^^ erg s“^ cm“^ with a photon index of 2 . 8 ~^ 2 t ^ 

luminosity < 3.0x10^® erg s“^. The SRGUT model implies a radio spectral index of 0.8j^Q^. 
There can be several explanations for this: (a) the energy going into particle acceleration is 
less, (b) the synchrotron losses and/or particle escape is more effective compared with the 
reverse shock yielding a lower X-ray flux for the forward shock zone, (c) the forward shock 
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is under-ionized as discussed in section (4.3) which results in a lower injection efficiency and 
thus less particles in the non-thermal regime. Ellison et ah (2004) shows that nonlinear 
diffusive shock acceleration in the forward and reverse shocks could evolve differently. The 
emission from the reverse shock may increase out to 1000 years and the efficiency diminishes 
after that together with the luminosity. The emission and efficiency from the forward shock 
is expected to increase slowly out to 40000 years and may be small for young SNRs (younger 
than 1000 yrs). 


7.2. The Classical Nova Remnants 

The studies of resolved old nova shells provide an opportunity of investigating several 
facts of the evolution of novae long after the initial eruption has subsided. Spectroscopy and 
imaging of the nebular remnants of old novae have been a diagnostic of physical conditions in 
(mode of excitation) and the morphology of the nova shells. Optical spectroscopy/imaging 
of nova shells has been widely used for this purpose. By direct imaging method, the ejecta 
of several old novae have been spatially resolved with evidence/detect ion of polar blobs and 
equatorial rings in many of them (Slavin et ah 1995; Gill & O’Brien 1998, 2000; Krautter et 
ah 2002; Bode 2002 and references therein). The expansion velocities of material in a nova 
shell are snggested to vary smoothly as a fnnction of position angle from a maximum in the 
polar regions to a minimum at the equator, which results in a prolate asymmetry (Lloyd 
et al. 1997). In addition to this, abundance gradients are expected as a result of TNRs 
(Thermo-Nuclear Runaway) in rotating oblate white dwarfs as detected from observations 
(RR Pic, DQ Her, V1370 Aql, V838 Her, etc.; Scott 2000). Recently, the HST imaging 
of HR Del indicates several emission knots in [OHI], [Nil] and Ha; particularly a bipolar 
structure in the first two compared with the latter (Harman & O’Brien 2003). Analysis of 
HR Del reveals that there are also indications of wind (from the central source) interaction 
with the old remnant shell affecting its evolution and structure. In general, optical and UV 
wavelength measurements of electron temperatnres in the nova shells indicate two groups of 
old remnants; hot (Te > 10*^ K) and cold (Tg < 10"^ K) nova shells (Williams 1992, Ferland 
1984). Some of the nova shells have been found to decelerate in time, such as DQ Her, GK 
Per, V603 Aql, T Pyx, and V476 Gyg, suggesting the existence of circnmstellar interaction 
(Duerbeck 1987). 

The first prediction of classical nova remnants as non-thermal radio sonrces has been 
made by Ghevalier (1977). However, there have not been any X-ray or radio wavelength 
(Bode et. al. 1987) detections of old classical nova remnants except for GK Per. The detailed 
theoretical framework on the hydrodynamical shock evolntion in classical nova remnants 
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is being developed (See Bode 2002 for a review). The circumstellar interaction has been 
modeled for recurrent novae (RS Oph: Bode & Kahn 1985; T Pyx: Contini & Prialnik 
1997). These models (esp. for T Pyx) predict the existence of a forward moving blast wave 
and a reverse shock moving into the ejecta as a result of colliding shells from the eruptions. 
Their model also predict a cold layer of material between the forward and reverse shocks 
(they do not study it in detail). The models indicate high electron densities of the order 
of a few hundred cm“^ in the shock zones. Another intriguing result from these studies is 
that the magnetic held is enhanced by factors of 10 to 100 due to instabilities and turbulent 
motion. 

The power from a simple shocked-shell model of thermal origin, as an X-ray emitting 
nebula, can be modeled assuming a 0.25 hlling factor for the volume of the shell (assumed 
to be spherical) and Ue ~ dUo (from strong shock conditions); 

Lt ~ 3.1 X (4) 

The temperature T is in units of 10^ K and the radius of the emitting shell R is in units 
of 3.1x10^® cm. The fast nova ejecta will start to cool as it sweeps up an equal amount of 
circumstellar material at Rcooz and tcooz- Once it sweeps up about 10 times its own mass, it 
will proceed into the Sedov phase; 

Rcooi ^ 0.03(M_5/no)^/^pc, (5) 

tcooi ^ 30(M_5/t;^ooono)^/^l/rs, (6) 

tsedov - 64(M_5/niooono)^'^^|/rs. (7) 

For Uo ~ 0.1-10 cm“^, kT ~ 0.1-10 keV and R Rcooi, the range of X-ray luminosity is a 
few X10^"^ to a few x 10^® ergs s“^ as the remnants start to cool. Assuming adiabatic remnants 
(at the onset of cooling), about 1% of this radiation should be detected. Thus, the maximum 
limit of X-ray emission is then about 10^^ ergs s“^ which is in good agreement with the X-ray 
luminosity of the remnant of GK Per (i.e., brightest remnant). The time to start to progress 
into the Sedov phase is < 96 yrs for GK Per, calculated using equation (7) which is consistent 
with the interpretation that the remnant is in the Sedov phase. I would like to note that the 
X-ray synchrotron emission from the remnant of GK Per is a special case and a result of the 
circumstellar environment and the magnetic helds. 

Theory would expect that the remnants should become fully radiative in the later stages 
of the shock evolution (i.e. Snowplow phase), increasing the luminosity in excess of xlO^® 
ergs s“^. However, by that time the remnant will expand to a much larger size where the 
decreasing surface brightness will affect this evolution and may no longer permit the remnant 
to be detected. 
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8. SUMMARY AND CONCLUSIONS 

The nova shell of GK Per is the hrst classical nova remnant resolved and detected 
in the X-ray wavelengths. Chandra Observatory has revealed the spectrum and detailed 
morphology of this remnant for the first time and this has helped to raise/resolve important 
questions related to classical nova eruptions, hydrodynamic flows and relativistic particle 
generation. In general, the nova remnant of GK Per evolves similar to young, type II SNRs 
and a shell is constructed via the interaction of the ejecta with the circumbinary medium 
which is a leftover wind and/or a Planetary Nebula. The X-ray nebula around GK per is 
brightest on the SW quadrant and toward West with a lumpy morphology. There is a central 
semi-circular region and ’’wing” like extensions toward SE and NW. The central circular 
region resembles the radio/VLA shell and the optical/HST images with several clumps/knots 
in the X-rays. Some features are collection of clumps and hlamentary structures, also, 
indicating radial elongations. 

The X-ray shell has a symmetric cone-shaped structure detected in a Ne IX emission 
line with a cone angle of about 74° to the central source which indicates the shock front. 
The shell is most likely expanding into a less dense medium in the NW and SE directions 
which comprises the ’’wings” with speeds (i.e., 2600 km s“^) more than twice the expansion 
speed in the SW direction (i.e., 1100 km s“^). The VLA/radio image of the remnant shows 
a trace along the ’’wings” detected in the X-rays as part of a larger diffuse radio emission 
stretching from NW to SE in the direction of the reflection nebulosity. The neon line emission 
and the neon wings detected toward the NW and SE, moving with speeds ~ 2600-2800 km 
s“\ originate from the clumped ejecta due to instabilities in the post-shock region. Such 
knots are expected to transit into the forward shock zone in time with possibly the aid of 
a circumstellar magnetic held and compression of the interaction region due to the effects 
of nonlinear diffusive shock acceleration (particle acceleration). It is also possible that some 
of these knots could have been created within the remnant or even at the onset of the 
nova explosion and have not signihcantly decelerated, eventually surpassing the shock zone 
particularly in the NW and SE directions where the interstellar density is low. The NE 
portion of the circumstellar environment does not show a distinct sign of interaction (i.e., a 
bright rim) as in the SW and lacks line emission. In general. X-ray nebula is the largest among 
the detected radio and optical shells. The shape of the X-ray nebulosity is a combination 
of the characteristics of the circumstellar medium around the nova and the geometrical 
projection effects that requires a special treatment with hydrodynamical models. Greation 
of a hemispherical X-ray shell, as detected for GK Per in this paper, is theoretically expected 
in regions with high ambient density gradients as revealed by studies in SNRs evolving in 
molecular clouds (e.g. W44, Gox et ah 1999). 



The nova shell of GK per is a site where energetic electrons are produced by diffusive 
shock acceleration (particle acceleration) conhrmed by the Chandra observation via the 
discovery of the synchrotron emission from the non-thermal electrons in the reverse shock 
zone. The detected power law emission in the X-rays reveals the same spectral index as in the 
radio wavelengths recovered with the fits using a SRCUT model of emission. The detected 
X-ray photon index is 2.3lo;9, spectral index is O.bSlgjg and the unabsorbed luminosity 
is 4.6 erg s“^. The remnant is age-limited. The non-thermal electron density is about 

3.01 xlO ^ (f=l) and the efficiency of injection in the reverse shock is about 1.4x10 This 
value is similar to the efficiency derived for SN 1006 which is the most efficient SNR in 
particle acceleration. However, the derived luminosity for GK Per is lower because the 
emitting volume is smaller. The cool X-ray temperature (0.1-0.3 keV) of the forward shock 
indicates that there is nonlinear diffusive shock acceleration occurring in this region, also. 
However, the efficiency and the luminosity is low (L2.<3xl0^® erg s“^ and r>2.8) and can 
not be resolved from the thermal emission. This complies with the theoretical expectations 
that for young SNRs the contribution from particle acceleration into the energetics of the 
systems will be small. 

The calculated break frequency for the remnant shows that GK Per is also a strong 
candidate for gamma-ray emission (out to the knee) due to the accelerated particles in the 
reverse shock region with E^ax about 15-30 TeV (this should be more or less correct for the 
forward shock region, too). 

The X-ray nebula has a low temperature component (below 2 keV), the hrst component, 
that is a (non-equilibrium ionization) thermal plasma close to ionization equilibrium, possibly 
dominated by emission lines, with kT~ 0.1-0.3 keV and an X-ray luminosity ~ 4.3x10^^ ergs 
s“^. This component shows the He-like Ne IX (i.e, most contribution may be coming from 
a forbidden emission line) and He-like N VI line emission derived from the spectral hts (and 
a possible O VH emission line). No H-like line emission is detected. This component is 
associated with the forward shock region and the transition zone owing to the low neutral 
Hydrogen column density. Thus, it is a mixed plasma of shocked circumstellar matter and 
ejecta material. The electron density within this region is 0.6-11.2 cm“^ (for a hlling factor 
of 1). The uq in the circumstellar medium is <60 cm“^ (no~4.1 cm“^). 

The higher energy component of the entire shell spectrum above 2 keV is most likely a 
plasma emission of bremsstrahlung origin with no resolved emission lines. The hts using a 
thermal bremsstrahlung model of emission reveal a shock temperature of kT = 1.041 q 5 keV 
and an X-ray luminosity (6.6-0.14) x 10^^ erg s“^. It is found to be heavily absorbed with 
Nh = (4.0-22.0) X10^^ cm“^ suggesting that there is a cold layer of material between the two 
X-ray components and the second component is most likely originating from the shocked gas 
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in the reverse shock zone. The X-ray analysis reveals that this cold layer is a plasma of solar 
abundances with N/Nq^II and Ne/Ne 0~18 where the high metalicity of the ejecta in this 
cooling layer is largely responsible for the high optical depth between 0.2 and 2 keV. It is 
also important to note that this high column density of neutrals absorbs the emission below 
1.5 keV where most of the line emission would be detected. 

The standard nova theory predicts enrichment of metal abundances produced prior to 
the outburst via mixing processes. This observation reveals for the first time a minimum 
enhancement of neon 13-21 and nitrogen 1-5 times its solar number fraction from a remnant 
in the X-rays that is 100 yrs old. Moreover, the existence of enhanced neutral abundances 
of neon and nitrogen in the ’’cold absorbing shell” supports this hnding. As a result, this 
conhrms that the outburst was a result of a violent TNR on an ONeMg WD, which would, 
also, explain the existence of a blast wave. The high nitrogen ratio is a direct evidence of 
the unstable nuclei which must have been ejected with the blast wave from the nova in 
the early stages. 

The kinetic energy density of the ejecta is (0.3-1.5) xl0“® erg cm“^ where the maximum 
limit comes from the calculated shocked mass (i.e. (2.1-38.5)x 10“^ M©)- The energy densi¬ 
ties of the thermal plasma, magnetic held and non-thermal electrons are 5.1x10“® erg cm“^, 
2.3x10“^® erg cm“^ and 3.0x10“^® erg cm“^, respectively, implying that the magnetic held 
and the non-thermal electrons are almost/in equipartition at the shock zone. The energy 
budget of the nova remnant indicates that not all the kinetic energy has appeared as ther¬ 
mal energy in the shock zone, in accordance with the fact that the remnant is young. Thus, 
considerable energy is still going into expansion of the shell and about less than 34% ap¬ 
pears as thermal energy with about less than 2% in the accelerated particles. The low X-ray 
luminosity (4.3x10^^ erg s“^) of the forward shock indicates that the remnant is adiabatic. 
The reverse shock is radiative. Overall, the nova shell of GK Per resemble a Sedov remnant 
(of SNRs). The cooling in the forward shock is also detected in the HST [Nil] images (very 
similar to Ha images) aligned with the X-ray data (see Figure 10a). 

The nova shell of GK Per remains to be unique with its characteristics over the entire 
electromagnetic spectrum. It behaves as a miniature SNR with typical characteristics, yet 
it is the hrst detection of GNRs (Glassical Nova Remnants), which can be a new class of 
X-ray emitters that have considerably low surface brightness/luminosity and shocked mass 
compared with SNRs; as a result evolving faster than, but similar to SNRs. The luminosity 
range for such objects are between a fewxlO^^ to a fewxlO^® ergs s“^ as the adiabatic 
remnants start to cool. After the outburst of a classical nova, the typical cooling timescale 
for the hard X-ray emission components due to wind-wind interactions is about three years 
(Balman et al. 1998). However, as the shell interacts with its circumstellar environment, it 



should produce about the range of luminosity above depending on how strong the shock will 
be (i.e., density contrast/compression ratio) and how the expansion will affect the surface 
brightness of the remnant. Two other candidates, DQ Her and RR Pic indicate residual 
emission coming from the location of their shells (Mukai, Still, & Ringwald 2003; Balman 
2002b). Recovering CNRs require simultaneous high sensitivity and spatial resolution and 
pose a challenge given the characteristics of the present and upcoming X-ray missions. The 
ACIS detector on Chandra and Epic PN on Newton-XMM could discover some of the close-by 
and luminous CNRs. 
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Table 1: Spectral Parameters of the Entire Spectrum of the Nova Shell for the Softer X-ray 
Component in the energy range 0.3-2 keV (ranges correspond to 2(7 errors) 


VPSHOCK §1 VMEKAL §2 BREMSS+3(GAUSSIAN) §3 


Njj (xlO^^ cm ^) 

U-oU_o.o7 

A 90 + 0.05 

n 10+018 
'-’•1^-0.16 

kT* (keV) 

n 1 1+0.04 
'4--L-L_0.02 

o.io+D? 

0’-+'-0.16 

Nitrogen (N/N©) 

0 1 + 1.6 

3.0t°;® 

N/A 

Oxygen (0/0©) 

1 8+1-'^ 

-^• 0 - 0.8 

N/A 

N/A 

Neon (Ne/Ne©) 

16-31^4 

9.23t^;^ 

N/A 

r (xlO^^ s cm“3) 

U.O_o .2 

N/A 

N/A 

EM (xl033 cm-3) 


66 7+133.3 
OO-' -63.24 

1 -1-238.4 

0+'4-51.9 

Line Energy (keV) 

N/A 

N/A 

Bremss: N/A 




Cl 

U.AZO_q Q 21 




G2 

A rc:7+0.005 

U.OOf _Q QQg 




G3 

0.907t°;r7 

Line Sigma (keV) 

N/A 

N/A 

Bremss: N/A 




G1 

0.005t°;°J° 




G2 

0.028t°;°?3 




G3 

0.003hO;OJ3 

Flux 


4 q+91 
u.»_ 2,2 

Bremss: 0. 32 ^ 0(03 




G1 

n 09 + 109.0 

U.OZ_o 28 




G2 

0.38tJ®2f 




G3 

0.032tjy5 

xl 

1.3 (84) 

1.5 (85) 

1.3 (87) 


§i§2§3rjn]^g spectral models are fitted simultaneously with the models labeled by the same number in Table 2. 

“1 keV - 1.2 xlO'^K. 

^Calculated using the normalization constant of the VMEKAL/VPSHOCK thermal plasma models (see 
sec[4.1]). A = (10“34y'47rD^) xEM where EM (Emission Measure) = f Ug nn dV (integration is over the 
emitting volume V). 

“Calculated using the normalization constant of the Thermal Bremsstrahlung model. 

A = (3.02x 10~35/47rD^) xE.M. where E.M. (Emission Measure) = f Ug nj dV (integration is over the 
emitting volume V). 

'^Unabsorbed Hard X-ray flux between 0.3 and 10 keV in units of xl0“^^ erg cm“^ s"^. The unconstrained 
parameters are kept at their best fit values when calculating model fluxes. 











Table 2: Spectral Parameters of The Entire Spectrum of the Nova Shell for the Harder 
X-ray Component above 2 keV (ranges correspond to 2a conhdence level errors) 


Model 

Nh 

(x 10^^ cm“^) 

kT, " 
(keV) 

T 

(x 10^^ s cm 

E.M. 

(xlO^'^ cm-3) 

Flux d 


NEpi 

PSHOCK§2 

BREMSS§3 

NEpi 

13.7l^-^ 

0.2t°;) s 

Q a+^-y 

0.96tJ'f 
1.04+J;f 
60.of 5]^ Q 

‘-'•'^^-0.089 

0.0-0.13^°-^^ f 
N/A 

6-7f4.3 

3.02ff5'^ 

9 9+11.1 
^■^-1.9 

9 4+6.4 

^■4-2.34 

0 001 +0 001 

U.UUi_Q QQQ5 

rtiy <^+322.8 
’^-25.8 

1 ’tD 

iOU.U_24g 6 

r <^+29.8 

^•^-4.45 

0.08t!];°6 

2.4 

Model 

Nh 

(x 10^^ cm“^) 

P 

a 

Break Erequency 
(xlO^® Hz) 

Flux d 

xl 

P0WER§3 

SRCUT§3t 

POWER^i 

10.511;^ 

0.29l° °4 g 

N/A 

0.04t°J 

N/A 

0.32+°-2® 

N/A 

N/A 

0 05+0 ® 

0’-UO_o.03 

N/A 

14000.0f^«™ 

6-2tf44 

0.3t4;°7 

i.94 


_ 00 

§i§2§3rp]^g spectral models are fitted simultaneously with the models labeled by the same number in Table 1. f 
The normalization of the SRCUT model (the radio flux at 1 GHz) is frozen at 0.024 Jy at the measured 
value (Reynolds & Chevalier 1984; Biermann et al. 1995) 

“1 keV ~ 1.2 xlO^RT. 

^Calculated using the normalization constant of the VMEKAL/PSHOCK model. A = (10“^^/47rD^)xE.M. 
where E.M. (Emission Measure) = f Ue nn dV (integration is over the emitting volume V). 

^Calculated using the normalization constant of the Thermal Bremsstrahlung model. 

A = (3.02x 10“^^/47rD^) xE.M. where E.M. (Emission Measure) = J Ue ni dV (integration is over the 
emitting volume V). 

'^Unabsorbed Hard X-ray flux between 0.3 and 10 keV in units of xlO“^^ erg cm“^ s“^. 

®The values for xt are the same as in Table 1 since the two emission models are fitted simultaneously. Thus, 
they will be omitted from this table. 

•^The range is derived using the values for the two parameters and tiow tiow is fixed at the value of 0.0 
for convenience, since it is the lower limit for the integration over the ionized volume and Tup shows the 
maximum limit on the ionization parameter at 2cr confidence level. 

®This value of Nh is derived from the fits with a single absorption parameter used for both X-ray spectral 
components. 



Table 3: Spectral Parameters of The Brightest Emission Region of the Nova Shell (ranges 


correspond to 2a conhdence level errors) 


Model for 


kTs =>• 

r 

Ne/Ne0 

N/Nq 

E.M. I" 

Flux 

xl 

Componentl 

(xlO^^ cm“^) 

(keV) 

(xlO^^ s cm ~^) 



(xlO®® cm“®) 



VPSHOCK§i 

n o't-|-U.U£) d 
-0.04 


0.0-2.8<2 4 ® 

y <:* + id.U 

1 Q + ^'^ 

' -0.2 

o.+lb 

f 

VMEKAL§2 

n oK+o.os d 
b'-^O_0 05 

r» -I +0.02 
^•^-0.01 

N/A 

7 0+7.7 
‘ •'^-4.1 

9 q+5.2 

1 1 y+15.0 

ii. / 2 

2 . 911 ^ 2 " 

f 

Model for 

Nh 

r 

a. 

Radio Flux 

Break Frq. 

Flux 

xl 


Component2 

(xlO^^ cm“^) 



at 1 GHz (Jy) 

(xlO^® Hz) 




POWERS^ 


9 Q-f-l.t) 

N/A 

N/A 

N/A 

1 y+ M.U 
-*-• ^-1.4 

1.15 


SRCUT§i 


N/A 

o.esini 

0.024 

1 9 + 2.3 
-*-•^-0.3 

1 ^7+125.3 
^-1.6 

1.16 


POWERS^ 

U.OO_o 23 

-0 2^^' ‘ 

N/A 

N/A 

N/A 

'5U.U_j^g 0 

1.9 


SRCUT§i 

^•^^-0.15 

N/A 

n oy+O.OO 

^•o<^_0.02 

0.024 

4.5< 

0.031“;°? 

3.0 


Model for 

Nh 

kTs 

r 

E.M. 

Flux 

xi 



Component2 

(xlO^^ cm“^) 

(keV) 

(xlO^^ s cm 

(xlO®® cm“®) 





PSHOCKS2 

A Q-f-H.U d 

20.9<i7,o 

96-95.9 

1 1 +14.b 

^•^-0.5 

0.8+iT < 

1.7 



MEKAL§1 

0.38tl\l g 

20 .0f J3 g 

N/A 

f) O+0-2 

U.O_o 2 


2.8 




§i§2The spectral models are fitted simultaneously with the models labeled by the same number. 

“1 keV ~ 1.2 xlO'^K. 

^Calculated using the normalization constant of the (V)MEKAL/(V)PSHOCK models. A = 
(10“^^/47rD^) xE.M. where E.M. (Emission Measure) = f Ug uh dV (integration is over the emitting volume 

V). 

“Unabsorbed Hard X-ray flux between 0.3 and 10 keV in units of xl0“^^ erg cm“^ s“^. The unconstrained 
parameters are kept at their best fit values when calculating model fluxes. 

'^This value of Nh is derived from the fits with two different absorption parameter used for the two X-ray 
spectral components. 

“The range is derived using the values for the two parameters and tiow is fixed at the value of 0.0 
for convenience, since it is the lower limit for the integration over the ionized volume and T„p shows the 
maximum limit on the ionization parameter at 2a confidence level. 

•^The values for xt s-re listed with the parameter results for Component 2 since the two emission models are 
fitted simultaneously. 

®This value of Njj is derived from the fits with a single absorption parameter used for both X-ray spectral 
components. 
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Fig. 1.— The X-ray nebula around the classical nova Persei (1901) between 0.3 and 10 
keV obtained with the Chandra ACIS-S (S3). The image shows emission 2a above the 
background. The resolution is 0".5 per pixel. It is smoothed using a variable Gaussian hlter 
with a=0".5-1". North is up and West is to the right. Central source is extracted from the 
image using a modeled PSF (see section 4.1) created with the aid of ChaRT, and MARX. 

Fig. 2.— The X-ray data of the shell htted with the VMEKAL-I-PSHOCK emission models 
including two different neutral hydrogen absorption. The smooth curve is the Chandra ACIS- 
S response to the composite model spectrum and the actual ACIS-S data are indicated with 
crosses. The dashed line shows the contribution of the VMEKAL model whereas the dot- 
dashed line indicates the contribution from the PSHOCK model of emission. The lower panel 
shows the residuals between the data and the model in standard deviations. In general, solar 
elemental abundances are assumed except for neon and nitrogen. 

Fig. 3.— (a) The source count rates extracted from an aperture of 12."5 versus the source 
count rates detected within an annular region from 12".5 to 67".5 in radius. It shows a scaling 
relation between the source count rates (largely piled-up) and the count rates derived from 
the PSF wings in the 2' vicinity of the pointed source. Both count rates are background 
subtracted, (b) The example calibration point source spectra used in Figure 3a. The central 
source count rate (pileup rate) increases from 0.19 to 0.21 c s“^ from row 1 column 1 to row 
2 col 2. The upper curves are the spectra from the source-region and the lower curves are 
spectra from the shell-region as denoted in the text. 

Fig. 4.— The top crosses (curve) indicate the data of the central source spectrum derived 
from an aperture of 12".5 (includes pileup photons). The middle crosses (curve) show the 
simulated central source spectrum scaled to match the expected count rate from the location 
covered by the X-ray shell which is 0.037 c s“^. The middle curve is the expected spectrum of 
the PSF wings for GK Per. The lower crosses (curve) is the un-piled central source spectrum 
(of GK Per) derived from the out-of-time events of the AGIS transfer streak. 

Fig. 5.— (a) The X-ray data of the shell cleaned from the effects of the PSF wings of 
the central source and htted with a single model of emission (VPSHOGK) together with 
a neutral hydrogen absorption model. The hgure stresses the signihcant deviation above 
1.6 keV. (b) The X-ray data of the shell cleaned from the effects of the PSF wings of the 
central source and htted with the VPSHOGK-I-NEI emission models including two diherent 
neutral hydrogen absorption models. The smooth curve is the Chandra AGIS-S response 
to the composite model spectrum and the actual AGIS-S data are indicated with crosses. 
The dashed line shows the contribution of the VPSHOGK model whereas the dot-dashed 
line indicates the contribution from the NEI model of emission. The lower panel shows the 
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residuals between the data and the model in standard deviations, (c) The same plot as in 
(b), but has only a single neutral hydrogen absorption parameter for both of the components. 
This spectrum is regrouped to have a Signal-to-Noise ratio of 13 per energy bin. 

Fig. 6.— (a) The figure shows several extraction regions used to derive the spectra to study 
the spectral variation within the remnant. The large circle is divided into 4 quadrants of SW, 
NW, NE, SE. The other three regions are constructed to study the brightest X-ray region. 
The extraction region is a PIE (annulus-sector as noted in the text) with (1) denoting the 
inner ridge, (2) the Peak X-ray zone, and (3) the outer ridge, (b) The X-ray data htted 
with the VMEKAL-I-PSHOCK emission models including neutral hydrogen absorption. The 
figure shows four spectra obtained by dividing the shell into four quadrants centered on the 
NE, SE, NW and SW with an inner radius of 12".5 and an outer radius of 67".5. 

Fig. 7.— (a) The X-ray data htted with a single model of emission (VPSHOCK) together 
with a single neutral hydrogen absorption model. This hgure stresses the signihcant deviation 
of the residuals above 1.6 keV and the existence of the second component, (b) The X-ray 
data htted with the VPSHOCK-I-SRCUT emission models including two diherent neutral 
hydrogen absorption. The smooth curve is the Chandra ACIS-S response to the composite 
model spectrum and the actual ACIS-S data are indicated with crosses. The dashed line 
shows the contribution of the VPSHOCK model whereas the dot-dashed line indicates the 
contribution from the SRCUT model of emission. The lower panel shows the residuals 
between the data and the model in standard deviations, (c) The same plot as in (b), but 
has only a single neutral hydrogen absorption parameter for both of the components. 

Fig. 8.— The image of the X-ray nebula extracted from 0.8 to 1.0 keV around the He-like 
Ne IX emission line. A suitable continuum between energy channels 1.12-1.4 (corresponding 
to the amount of energy channels within the line at the continuum level) is subtracted from 
the image. The resolution of the image is 1" per pixel. North is up and West is to the 
right. The image is smoothed using a Gaussian hlter with a=2". The cross point shows the 
location of the central source. 

Fig. 9.— (a) The left-hand panel is the image of the X-ray nebula in the 0.25-0.50 keV 
energy band, (b) The middle panel is the image of the X-ray nebula in the 0.51-1.50 keV 
energy band. The two images in (a) & (b) show two sigma emission above the background 
at l"/pixel resolution, smoothed using a Gaussian hlter with a variable a=V'-2". (c) The 
right-hand panel is the image of the X-ray nebula in the 1.51-8.0 keV energy band at a lower 
resolution of 3"/pixel (compared with (a) & (b)) smoothed using a Gaussian hlter with a 
(7=6". The image in (c) shows three sigma emission above the background. The overlay are 
contours of the 0.51-1.50 keV image produced at the same spatial resolution and smoothed 



using the same Gaussian. 


Fig. 10.— (a) The figure shows the intensity contours of the HST image (1997) obtained with 
the [Nil] filter overlaid on the X-ray image of shell of GK Per. Both the HST contours and 
the X-ray image show emission 2a above the background. The resolution of the combined 
image is 2" per pixel. Both images are also smoothed by a Gaussian filter with a=2". 
The [Nil] intensity contours are based on histogram eqnalization starting with 0.6a to 68(T 
increments between contonr levels, (b) The figure shows the radio intensity contours at 
1.425 GHz overlaid on the X-ray image of the shell of nova Persei (1901). The X-ray image 
and the radio contours show emission 2a and 1.5a above the backgronnd, respectively. The 
resolution of the combined image is 2" per pixel. Both images are also smoothed by a 
Gaussian filter with a=\". The radio intensity contours are constructed with linear 1.4a 
increments between contonr levels. 



03^31'^20'' 03'"31'^15® 03^3l'^10'' 03^3l'^05 


41 



2.2x10 ^ 2.6x10 ^ 3.5x10 ^ 5.7x10 ^ 6x10 





42 



(keV) 















43 



(s/o) snjpEJ u! |0X!d gsi. oi 93 ujojj uo!60j jeinuue ue u\ ^ve^ lunoQ 


0.01 












/\8>l/o8s/s:|.unoo pazipLUJOu 












Declination 



Right ascension 






- 50 - 



- (N 


in 

d 


channel energy (keV) 






























































































































































































































































3:J1:20 


15 


10 


05 


Right ascension 




(0.36-0.6 keV) (0.61-1.6 keV) (1.61-8.0 keV) 


55 












56 











